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A  SIMPLE  MODEL  FOR  DISSOCIATION  OF  DIATOMIC  MOLECULES 


1  Equation  of  state 

The  equation  of  state  (EOS)  which  is  currently  used  in  the  NRL  hydro  code  is  an 
improved  version  of  the  EOS  described  in  reference  [1].  In  this  section  We  shall 
briefly  review  the  parts  of  this  EOS  which  are  relevant  to  dissociation. 

1.0.1  The  basic  approach 

The  fundamental  assumption  in  the  formulation  of  this  EOS  is  that  the  free 
energy  of  the  system  (and  therefore  all  the  thermodynamic  quantities)  may  be 
written  as  a  superposition  of  three  terms  representing  the  free  energy  associated 
with; 

1.  atomic  and  electronic  interactions  at  zero  tCTQperatur6-E;(p) 

2.  Thermal  motion  of  atoms  and  ions-Fn(p,  T) 

3.  Thermal  motion,  excitation  and  ionization  of  electrons.-Fe(p.  T) 

(  r  is  the  temperature  and  p  is  the  mass  density.) 

1.0.2  Nuclear  contribution  to  the  EOS 

In  order  to  describe  the  motion  of  atoms  and  ions  in  solid,  liquid  and  gaseous 
states  in  a  sin^e  analytical  formula,  the  authors  of  [1]  use  an  interpolation 
method.  Equations  (4.18), (4.19), (420)  of  reference  [1]  present  an  interpolation 
between  the  free  energy  of  solid  according  to  the  Debye-Griineisen  model  and  the 
free  energy  for  ideal  gas.  Elxplicdtly  these  equations  read; 


Fn  =  NoksT |31n(l  - e"^)  - 1  +  |ln(l  +  ^)| 

(1) 

where 

(2) 

2‘KkB  0^ 

In  these  equations,  h  is  the  Plank  constant,  ks  is  the  Boltzman  constant,  Ni  is  the 
number  of  atoms  per  unit  mass  with  atomic  number  Zi  and  atomic  mass  mj,  Nq 
is  the  number  of  atoms  per  unit  mass;  No  =  Y^Ni.  0  is  the  Debye  temperature 

which  in  this  model  has  the  following  functional  dependence  on  the  mass  density: 
(equation  (4.12)  in  reference  [1]): 
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In  the  above  formula,  Oq  and  To  are  reference  values  of  the  Debye  temperature 
and  the  Griineizen  coefficient  at  a  reference  mass  density  po  .The  pr^sme  and 

energy  associated  with  this  free  energy  (by  Pn  —  ~  sr  ) 


are: 


Pn  =  pNoksT 


3r  +  ^ 

l  +  'J' 


(4) 


and 


En  =  ^NokeT 


2  +  $ 
1  +  ^ 


where  F  is  the  Griineisen  coefficient: 


r  =  i.— 

6  dp 


(5) 


2  The  partition  function 

In  the  present  section  We  will  derive  the  effective  ’’single  particle  partition  func¬ 
tion”  from  which  the  free  energy  of  equation  1  ,and  therefore  also  the  nuclear 
parts  of  the  equation  of  state  used  at  NRL  are  derivable.  In  section  3  We  will 
use  this  partition  function  in  the  derivation  of  the  equation  for  the  degree  of 
dissociation. 

Consider  the  function  ^  of  equation  2 


Nlh^p^Tyrl 

2'irkB  N^rnf, 

Using  the  relation  Nq  =  5^  iVj  in  the  last  equation,  one  gets 


(6) 


(7) 
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Denote  the  number  of  atoms  with  atomic  number  Zj  and  atomic  mass  m/  in 
the  volume  V  by  rji  and  the  total  number  of  atoms  by  r\  (i.e.  rj  =  P  = 
i  X^T?imi,and  Ni  =  =  F?)  •  Using  this  notation  and  rearranging  formula 

7  ,  '9  may  be  rewritten  as: 


Tlie  tom  in  the  denominator  of  the  right  hand  side  of  equation  8  is  the  well 
known  single  particle  partition  function  for  noninteracting  particles  of  mass  mi 
in  a  volume  V  (ideal  gas) 


qr=v 


(9) 


The  term  in  the  numerator  is  the  partition  function  for  solid  (Debye  Griineisen 
model,  in  the  limit  |:  «  1) 


n 


Return  now  to  equation  1,  in  the  limit  |  «  1,  it  may  be  written  as: 
F„=NotBr  jln(i(|;)?) +|h.(l  +  4’)} 


2i.> 


^NokeT 


+  ln  1  + 


n(*? 


9?/ 


=win|(n(^)")[i+n(!j)*’]“} 
[[i+n(»f)"] 


=  NokBT\ii{ 


i 


(10) 

(11) 


(12) 
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=  -NokBT\n{ 


ZUL 


n  («”)’■' 

I 


We  therefore  conclude  that  the  interpolated  free  energy  of  equation  1  is  derivable 
from  the  partition  function: 


In  the  limit  ^  »  1,  the  interpolated  partition  function  reduces  to  the  partition 
function  for  a  mixture  of  ideal  gases; 


In  the  opposite  limit.'®' 


^-►OO  -*■*  \t]l  J 

0  the  partition  function  reduces  to  the  Debye  formula: 


limZ„  =  n(9'')" 


p-^OO 


(Since  X)  = »?) 

NotiU  that  far  from  these  limits,  the  partition  function  13  does  not  factor 
into  the  partition  functions  the  various  species  participating  in  the  mixture.. 

3  Dissociation  of  diatomic  molecules 


3.1  Law  of  mass  action 

The  incorporation  of  dissociation  process  into  the  equation  of  state  is  based  on 
the  law  of  mass  action[2].  According  to  this  law,  the  equihbrium  constant  for 
dissociation  of  diatomic  molecules  2X)  is  related  to  the  single  particle 

partition  function  for  atoms,  mid  for  molecules,  qx^i  by. 


•qxJV  qxjy 


(14) 


Where  rjx  is  the  number  of  atoms  in  the  volume  V,  t]x^,  is  the  number  of 
molecules  in  the  volume  V.  Equation  14  together  with  the  obvious  relation 


f]x  "b  2^*2  ^ 


(15) 
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(which  states  that  the  total  number  of  atoms  in  the  system  is  constant)  allows  us 
to  solve  for  the  densities  of  the  atoms  and  the  molecules  jt®  =  Vxl^  and  n®2  — 
in  terms  of  the  partition  functions,  qJV  and  q^jV.  Strictly  speaking,  the 
law  of  mass  action  is  applicable  only  to  mixtures  of  reacting  ideal  gases,  in  the 
present  work  We  will  ignore  this  restriction  (without  any  justification)  and  use 
equation  14  together  with  the  interpolated  partition  function. 

3.2  Single  particle  interpolated  partition  function  for  atoms 
and  molecules 

By  equation  13, the  single  particle  partition  function  is: 


(16) 


This  partition  function  takes  into  account  only  translation  of  atoms  and  mole¬ 
cules.  For  molecules,  rotation ,  vibration  and  bonding  energy  of  the  molecule  also 
contribute.  In  the  next  section  We  shall  use  qn  as  the  smgle  particle  partition 
function  for  atoms,  for  molecules  We  shall  multiply  by  the  partition  functions  for 
vibration  and  rotation  and  dissociation,  q^  ,q^ » i*®-j  hi  equation  14  We  shall 


use: 


and; 


g®2  =  qn{2m)(^  (^qdiss  = 


1  + 


— wQ^q^Qdiss 
2 


(18) 


Strictly  speaking,  the  right  hand  side  of  eqiration  17  should  be  multiplied  by  the 
partition  function  for  electrons  in  an  atom  and  the  right  hand  side  of  equation 
18  hy  the  partition  function  of  electrons  in  a  molecule.  In  the  present  work  a  less 
rigorous  approach  is  used,  the  effect  of  electrons  is  accounted  for  only  trough  the 
dissociation  partition  function; 


QdisB  —  ^ 


(19) 
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which  represent  the  ratio  between  the  partition  function  of  two  dissodated  atoms 
and  that  of  a  molecule  in  a  phenomenological  way  with  a  density  dependent 
dissociation  energy.  We  use  the  function  [3]: 

fo  v>vA 

E^u.  =  E^..o  •  I  V<V^j 

For  We  shall  take  the  partition  function  of  harmonic  oscillator  in  the  limit 
where  the  temperature  is  higher  than  the  vibrational  temperature  By  =  ^  « 
T; 


n=0 


==e  2t 


-en/ksT  _ 

n=d) 

oo 

\  e  T  ^  e 

n=0 

1 

Ll  —  e  r  J 

(20) 


2sinh(|^) 

For  q^  we  shall  take  the  partition  function  for  rotator  with  moment  of  inertia 


q^  = 

/=o 


y^(2Z  +  l)e  ~  /*  C  ®  ^ 

/=o  0 


Sir^IkeT 


where: 


_r 

(21) 


$R  = 


2IkB 


(22) 


(For  hydrogen  molecules  the  vibrational  temperature  is  9v  =  6320  fC,  the  rota¬ 
tional  temperature  is  9jt  —  87.5  K,  and  dissociation  energy  is  Ejnss  4.48  eV^.) 


3.3  An  equation  for  the  densities  of  molecules  and  atoms 

Consider  a  mixture  of  tj,  atoms  of  mass  m,  t?®,  diatomic  molecules  of  mass  2m 
and  rii  inert  atoms  of  mass  mi  in  a  volume  V.  At  every  given  temperature  T  and 
mass  density  p,  that  the  dissociation  reaction  has  arrived  a  steady  state, 

the  number  of  atoms  and  molecules  in  every  volume  element  is  determined  by  the 
law  of  mass  action  14.  In  the  present  section  We  shall  use  the  law  of  mass  action 
14  and  the  interpolated  partition  functions  derived  in  section  3.2  and  derive  an 
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Equations  30  and  27  are  two  equations  for  the  atomic  density  and  the 
molecular  density  n*2  for  given  values  of  mass  density  p  and  temperature  T. 

In  the  rest  of  the  section  We  shall  simplify  the  notation  of  these  equations. 
Define  the  quantities; 


I3{p,T)=2 


q^{T)q^iT) 

1  (2'KmkftT\2 

X  V  h®  / 


and 


Tlx 


In  terms  of  these  quantities,  equation  30  takes  the  form: 


(32) 

(33) 


8 


(34) 


+(¥)’]  . 

Ri[i  +  Rif 

3.4  The  procedure  for  obtaining  the  atomic  and  molecu¬ 
lar  densities  for  a  given  state  with  known  total  mass 
density  p  and  temperature  T. 

For  given  temperature  T  and  mass  density  p ,  the  functions  ^(p)>  ^(p)>  C(p>  ^)>  0 

may  be  evaluated  e35)li(atly  by  using  formulas  3,26  33,  31.  Thus,  the  only  un¬ 
known  in  equation  34  is  R.  Solving  equation  34,  one  gets  i2(p,  T)  ,  the  atonuc 
and  molecular  densities  Ux  and  71*3  are  obtained  by: 

Tlx  =  RX  (35) 

(36) 

In  section  4  We  will  describe  a  numerical  analysis  of  the  equations  and  an  example 
where  the  above  procedme  is  applied  to  the  case  where  ^  «  1.  In  this  case  the 
equations  may  be  solved  analyticahy.  In  section  5  a  Fortran  routine  for  applying 
the  procedure  and  solving  for  the  atomic  and  molecular  densities  in  the  general 
case  is  presented. 


4  Numerical  Examples 

4.0.1  constants  and  functions  definition 


h  =  erg-sec 

fcfi  =  1.38  ♦  10““  erg/kelvin 
for  Deuteriiun 

m  =  2  ♦  1.67  *  10-24,  m  =  3. 34  X  lO'^^ 

Or  =  87.5  k 
Ov  =  6320  k 

EdUa  =  4.48  eV=  4.48  *  1.6022  *  10““  :  erg  i.e.  =  7. 1779  x  10  erg 
take: 


Fo  =  ffo  =  1*8 
do  =  240  Kelvin 
po  =  0.3017  g/cc 
mi  =  16  *  1.67  *  10“24 
mi  =2. 672x10-23 


a  =  0. 
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2sinh\^) 

MT)  =  4 


/3(p,r)  =  2 


gv(T)qv(T) 


aP,T) 


(tn+a*miy 
/2'trmkfiT\ 

!=  — .  »■  ^  n 


(  2'frmkj^T  \  ] 

{—far-) 


2  ^difia 
g3  kffT 


4.0.2  Debye  temperature 

Let  us  start  by  plotting  the  formula  for  Debye  temperature  used  in  reference  [1] 


4.0.3  The  interpolation  factor,  ^ 

Notice  that  at  low  densities,  0  becomes  very  low,  of  the  order  of  few  degrees. 
This  in  turns  means  that  ^  becomes  very  small.  For  a  single  type  of  atoms, 

^  =  (^V  ,  when  this  ratio  becomes  much  smaller  than  1,  as  can  be  seen 
imme^S;dy  from  equation  13,  the  partition  function  and  consequently  ^  the 
other  thermodynamics  quantities,  reach  the  limit  of  ideal  gas.  The  surface  4  (p,  T) 
is  shown  in  the  figure  below.  The  line,  in  the  p,  T  plane,  formed  by  the  intersection 
between  the  ^  surface  and  the  bottom  of  the  box  =  0.2)  determines  the  region 
for  which  ^  «  1  and  the  ideal  gas  behavior  prevails.  Similarly  ,  -  ,  , 

The  line  formed  by  the  intersection  between  the  ^  surface  and  the  top  of  the  box 
(e  =  5)  determines  the  region  for  which  ^  »  1  and  soUd  state  behavior  prevaik. 
In  the  region  between  these  lines  the  full  interpolation  formulas  1,4  and  5  sho 
be  used. 
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The  interpulation  factor  ^  as  a  function  of  p{gfcc)^T{Kelvin) 


4.0.4  Ideal  gas  case 

In  the  following  We  shall  consider  the  case  ^  «  1.  This  simple  analytically 
solvable  case  will  serve  as  a  test  case  when  solving  numerically  for  the  general 
case.  When  ^  «  1,  equation  34  takes  the  form: 


or, 


R^  +  -^R- 

2/?f 


The  solutions  for  this  equation  is: 


(37) 


(38) 


R(p,T) - i 


(39) 


Only  the  upper  sign  solution  is  physical  since  this  is  the  only  solution  which 
leads  to  values  in  the  correct  range,  0  <  R  <  1  .  The  figure  below  shows  the 
dissociation  ratio  as  a  function  of  p  and  T. 
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and  at  higher  temperatures: 


At  the  high  density  comer  ot  the  box  «0.13, 3000)  =  .  26644  so  beyond  this 
density  the  ideal  gas  approximation  is  not  vaUd  anymore  and  we  have  to  solve 
the  more  general  equation  34 

5  The  general  case:  A  routine  for  solving  for  atomic  and 
molecular  densities 

5.1  The  program 

The  program  for  solving  the  equation  for  the  molecular  density  and  the  atomc 
density  n,,  for  the  general  case,  (equation  34)  is  presented  in  the  Appendix.  For 
given  values  of  a,  p  and  T  (a,rho  and  t  in  the  program),  the  program  e^^uat^ 
the  effective  density  of  hydrogen  atoms  A(p)  and  the  coefficients  ^(p,r)  and 
l3(p,T)  (csi  and  beta  in  the  program)  of  equation  34  .  Having  these  coefficien  s 
the  program  calls  the  routine  rtsafe  which  solves  equation  34  for  R.  (The  routine 
rtsafe  is  based  on  the  Newton-Raphson  method  for  finding  zeros  of  functions). 
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the  molecular  density  and  the  atomic  density  71*2  evaluated  by  using 
equations  35  and  36. 

5.2  An  example 

The  routine  described  above  was  incorporated  into  the  EOS.  Figure  (5)  describes 
the  results  of  the  application  of  the  improved  EOS  to  measurements  described 
in  reference  [4]  .  The  squares  describe  the  experimental  Hugoniot.'pie  left  line 
(”ion”)  describes  the  Hugoniot  curve  with  the  EOS  without  dissociation.  The 
next  Hugoniot  (”sim”)  was  obtained  by  using  EOS  with  the  assumption  that 
dissociation  fraction  is  just  ^{-Edu^lksT)  .  The  curve  ”mol”  was  obtained 
by  using  EOS  with  the  model  described  above  but  with  a  constant  dissociation 
energy.  The  "new  model”  result  was  obtained  by  using  the  EOS  with  the  complete 
model.  The  line  "ref  11”  corresponds  to  the  model  of  reference  (11)  in  [4]. 


0 

1.. 

3 
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Figure  (5):  The  results  of  the  application  of  the  improved  EOS 
to  the  measurements  described  in  reference  [4]. 


13 


References 

[1]  S.  L.  Thompson  and  H.  S.  Lauson  Improvements  in  the  chart  D  Radiation- 
hydrodynamic  code  ///.'  revised  analytic  equation  of  state.  Report  sc-RR-71 

0714  (1975).  . 

[2]  Chang  L.  Tien  and  John  H.  Lienhard  Statistical  Thermodynamics  (McGraw- 

Hill)  ,  ^ 

[3]  David  A.  Young  and  Ellen  M.  Corey  J.  Appl.  Phys  78,  3748  (1995). 

[4]  L.  B.  Da  Silva  et.  al.  Phys.  Rev.  Letters  78,  483  (1997). 


14 


